SLAC-PUB-13556 



Search for the Rare Leptonic Decays B + — > l + V(_ (£ = e, fx) 

B. Aubert, 1 Y. Karyotakis, 1 J. P. Lees, 1 V. Poireau, 1 E. Prencipe, 1 X. Prudent, 1 V. Tisserand, 1 J. Garra Tico, 2 
E. Grauges, 2 M. Martinelli ab , 3 A. Palano ab , 3 M. Pappagallo ab , 3 G. Eigen, 4 B. Stugu, 4 L. Sun, 4 M. Battaglia, 5 

D. N. Brown, 5 L. T. Kerth, 5 Yu. G. Kolomensky, 5 G. Lynch, 5 I. L. Osipenkov, 5 K. Tackmann, 5 T. Tanabe, 5 
C. M. Hawkes, 6 N. Soni, 6 A. T. Watson, 6 H. Koch, 7 T. Schroeder, 7 D. J. Asgeirsson, 8 B. G. Fulsom, 8 C. Hearty/ 

T. S. Mattison, 8 J. A. McKenna, 8 M. Barrett, 9 A. Khan, 9 A. Randle-Conde, 9 V. E. Blinov, 10 A. D. Bukin, 10 ^Q 
A. R. Buzykaev, 10 V. P. Druzhinin, 10 V. B. Golubev, 10 A. P. Onuchin, 10 S. I. Serednyakov, 10 Yu. I. Skovpen, 10 

E. P. Solodov, 10 K. Yu. Todyshev, 10 M. Bondioli, 11 S. Curry, 11 I. Eschrich, 11 D. Kirkby, 11 A. J. Lankford, 11 
P. Lund, 11 M. Mandelkern, 11 E. C. Martin, 11 D. P. Stoker, 11 S. Abachi, 12 C. Buchanan, 12 H. Atmacan, 13 

J. W. Gary, 13 F. Liu, 13 O. Long, 13 G. M. Vitug, 13 Z. Yasin, 13 L. Zhang, 13 V. Sharma, 14 C. Campagnari, 15 
T. M. Hong, 15 D. Kovalskyi, 15 M. A. Mazur, 15 J. D. Richman, 15 T. W. Beck, 16 A. M. Eisner, 16 C. A. Heusch, 16 
J. Kroseberg, 16 W. S. Lockman, 16 A. J. Martinez, 16 T. Schalk, 16 B. A. Schumm, 16 A. Seiden, 16 L. Wang, 16 
L. O. Winstrom, 16 C. H. Cheng, 17 D. A. Doll, 17 B. Echenard, 17 F. Fang, 17 D. G. Hitlin, 17 I. Narsky, 17 
T. Piatenko, 17 F. C. Porter, 17 R. Andreassen, 18 G. Mancinelli, 18 B. T. Meadows, 18 K. Mishra, 18 M. D. Sokoloff, 18 

P. C. Bloom, 19 W. T. Ford, 19 A. Gaz 19 J. F. Hirschauer, 19 M. Nagel, 19 U. Nauenberg, 19 J. G. Smith, 19 
S. R. Wagner, 19 R. Ayad, 20 -[l| A. Soffer, 20 '[! W. H. Toki, 20 R. J. Wilson, 20 E. Feltresi, 21 A. Hauke, 21 H. Jasper, 21 
T. M. Karbach, 21 J. Merkel, 21 A. Petzold, 21 B. Spaan, 21 K. Wacker, 21 M. J. Kobel, 22 R. Nogowski, 22 
K. R. Schubert, 22 R. Schwierz, 22 A. Volk, 22 D. Bernard, 23 G. R. Bonneaud, 23 E. Latour, 23 M. Verderi, 23 
P. J. Clark, 24 S. Playfer, 24 J. E. Watson, 24 M. Andreotti Qb , 25 D. Bettoni a , 25 C. Bozzi a , 25 R. Calabrese ab , 25 
A. Cecchi ab , 25 G. Cibinetto ab , 25 E. Fioravanti afc , 25 P. Franchini ab , 25 E. Luppi afc , 25 M. Munerato afc , 25 M. Negrini ab , 25 

A. Petrella ab , 25 L. Piemontese a , 25 V. Santoro ab , 25 R. Baldini-Ferroli, 26 A. Calcaterra, 26 R. de Sangro, 26 

G. Finocchiaro, 26 S. Pacetti, 26 P. Patteri, 26 I. M. Peruzzi, 26 ! M. Piccolo, 26 M. Rama, 26 A. Zallo, 26 R. Contri Qb , 27 
E. Guido, 27 M. Lo Vetere ab , 27 M. R. Monge ab , 27 S. Passaggio a , 27 C. Patrignani ab , 27 E. Robutti Q , 27 S. Tosi ab , 27 
K. S. Chaisanguanthum, 28 M. Morii, 28 A. Adametz, 29 J. Marks, 29 S. Schenk, 29 U. Uwer, 29 F. U. Bernlochner, 30 

V. Klose, 30 H. M. Lacker, 30 D. J. Bard, 31 P. D. Dauncey, 31 M. Tibbetts, 31 P. K. Behera, 32 M. J. Charles, 32 
U. Mallik, 32 J. Cochran, 33 H. B. Crawley, 33 L. Dong, 33 V. Eyges, 33 W. T. Meyer, 33 S. Prell, 33 E. I. Rosenberg, 33 
A. E. Rubin, 33 Y. Y. Gao, 34 A. V. Gritsan, 34 Z. J. Guo, 34 N. Arnaud, 35 J. Bequilleux, 35 A. D'Orazio, 35 
M. Davier, 35 D. Derkach, 35 J. Firmino da Costa, 35 G. Grosdidier, 35 F. Le Diberder, 35 V. Lepeltier, 35 A. M. Lutz, 35 

B. Malaescu, 35 S. Pruvot, 35 P. Roudeau, 35 M. H. Schune, 35 J. Serrano, 35 V. Sordini, 35 '0 A. Stocchi, 35 
G. Wormser, 35 D. J. Lange, 36 D. M. Wright, 36 I. Bingham, 37 J. P. Burke, 37 C. A. Chavez, 37 J. R. Fry, 37 

E. Gabathuler, 37 R. Garnet, 37 D. E. Hutchcroft, 37 D. J. Payne, 37 C. Touramanis, 37 A. J. Bevan, 38 C. K. Clarke, 38 

F. Di Lodovico, 38 R. Sacco, 38 M. Sigamani, 38 G. Cowan, 39 S. Paramesvaran, 39 A. C. Wren, 39 D. N. Brown, 40 

C. L. Davis, 40 A. G. Denig, 41 M. Fritsch, 41 W. Gradl, 41 A. Hafner, 41 K. E. Alwyn, 42 D. Bailey, 42 R. J. Barlow, 42 

G. Jackson, 42 G. D. Lafferty, 42 T. J. West, 42 J. I. Yi, 42 J. Anderson, 43 C. Chen, 43 A. Jawahery, 43 

D. A. Roberts, 43 G. Simi, 43 J. M. Tuggle, 43 C. Dallapiccola, 44 E. Salvati, 44 S. Saremi, 44 R. Cowan, 45 D. Dujmic, 45 
P. H. Fisher, 45 S. W. Henderson, 45 G. Sciolla, 45 M. Spitznagel, 45 R. K. Yamamoto, 45 M. Zhao, 45 P. M. Patel, 46 

S. H. Robertson, 46 M. Schram, 46 A. Lazzaro ab , 47 V. Lombardo a , 47 F. Palombo ab , 47 S. Stracka ab , 47 
J. M. Bauer, 48 L. Cremaldi, 48 R. Godang , 48 'B R- Kroeger, 48 P. Sonnek, 48 D. J. Summers, 48 H. W. Zhao, 48 
M. Simard, 49 P. Taras, 49 H. Nicholson, 50 G. De Nardo ab , 51 L. Lista a , 51 D. Monorchio ab , 51 G. Onorato ab , 51 
C. Sciacca ab , 51 G. Raven, 52 H. L. Snoek, 52 C. P. Jessop, 53 K. J. Knoepfel, 53 J. M. LoSecco, 53 W. F. Wang, 53 
L. A. Corwin, 54 K. Honscheid, 54 H. Kagan, 54 R. Kass, 54 J. P. Morris, 54 A. M. Rahimi, 54 J. J. Regensburger, 54 
S. J. Sekula, 54 Q. K. Wong, 54 N. L. Blount, 55 J. Brau, 55 R. Frey, 55 O. Igonkina, 55 J. A. Kolb, 55 M. Lu, 55 
R. Rahmat, 55 N. B. Sinev, 55 D. Strom, 55 J. Strube, 55 E. Torrence, 55 G. Castelli ab , 56 N. Gagliardi ab , 56 
M. Margoni ab , 56 M. Morandin a , 56 M. Posocco a , 56 M. Rotondo a , 56 F. Simonetto ab , 56 R. Stroili ab , 56 C. Voci ab , 56 
P. del Amo Sanchez, 57 E. Ben-Haim, 57 H. Briand, 57 J. Chauveau, 57 O. Hamon, 57 Ph. Leruste, 57 G. Marchiori, 57 
J. Ocariz, 57 A. Perez, 57 J. Prendki, 57 S. Sitt, 57 L. Gladney, 58 M. Biasini Qb , 59 E. Manoni ab , 59 C. Angelini ab , 60 

G. Batignani ab , 60 S. Bettarini ab , 60 G. Calderini ab , 60 '[l!] M. Carpinelli ab , 60 '[il A. Cervelli ab , 60 F. Forti ab , 60 
M. A. Giorgi ab , 60 A. Lusiani ac , 60 M. Morganti ab , 60 N. Neri ab , 60 E. Paoloni ab , 60 G. Rizzo ab , 60 J. J. Walsh a , 60 



2 



D. Lopes Pegna, 61 C. Lu, 61 J. Olsen, 61 A. J. S. Smith, 61 A. V. Telnov, 61 F. Anulli a , 62 E. Baracchini afc , 62 
G. Cavoto Q , 62 R. Faccini Qb , 62 F. Ferrarotto , 62 F. Ferroni afc , 62 M. Gaspero ab , 62 P. D. Jackson a , 62 L. Li Gioi a , 62 
M. A. Mazzoni a , 62 S. Morganti a , 62 G. Piredda a , 62 F. Renga ab , 62 C. Voena a , 62 M. Ebert, 63 T. Hartmann, 63 
H. Schroder, 63 R. Waldi, 63 T. Adye, 64 B. Franek, 64 E. O. Olaiya, 64 F. F. Wilson, 64 S. Emery, 65 L. Esteve, 65 

G. Hamel de Monchenault, 65 W. Kozanecki, 65 G. Vasseur, 65 Ch. Yeche, 65 M. Zito, 65 M. T. Allen, 66 D. Aston, 66 
R. Bartoldus, 66 J. F. Benitez, 66 R. Cenci, 66 J. P. Coleman, 66 M. R. Convery, 66 J. C. Dingfelder, 66 J. Dorfan, 66 

G. P. Dubois-Felsmann, 66 W. Dunwoodie, 66 R. C. Field, 66 A. M. Gabareen, 66 M. T. Graham, 66 P. Grenier, 66 

C. Hast, 66 W. R. Innes, 66 J. Kaminski, 66 M. H. Kelsey, 66 H. Kim, 66 P. Kim, 66 M. L. Kocian, 66 D. W. G. S. Leith 66 
S. Li, 66 B. Lindquist, 66 S. Luitz, 66 V. Luth, 66 H. L. Lynch, 66 D. B. MacFarlane, 66 H. Marsiske, 66 R. Messner, 66 'Q 
D. R. Muller, 66 H. Neal, 66 S. Nelson, 66 C. P. O'Grady, 66 I. Ofte, 66 M. Perl, 66 B. N. RatclifT, 66 A. Roodman, 66 
A. A. Salnikov, 66 R. H. Schindler, 66 J. Schwiening, 66 A. Snyder, 66 D. Su, 66 M. K. Sullivan, 66 K. Suzuki, 66 
S. K. Swain, 66 J. M. Thompson, 66 J. Va'vra, 66 A. P. Wagner, 66 M. Weaver, 66 C. A. West, 66 W. J. Wisniewski, 66 
M. Wittgen, 66 D. H. Wright, 66 H. W. Wulsin, 66 A. K. Yarritu, 66 K. Yi, 66 C. C. Young, 66 V. Ziegler, 66 X. R. Chen, 67 

H. Liu, 67 W. Park, 67 M. V. Purohit, 67 R. M. White, 67 J. R. Wilson, 67 P. R. Burchat, 68 A. J. Edwards, 68 
T. S. Miyashita, 68 S. Ahmed, 69 M. S. Alam, 69 J. A. Ernst, 69 B. Pan, 69 M. A. Saeed, 69 S. B. Zain, 69 S. M. Spanier, 70 
B. J. Wogsland, 70 R. Eckmann, 71 J. L. Ritchie, 71 A. M. Ruland, 71 C. J. Schilling, 71 R. F. Schwitters, 71 
B. C. Wray, 71 B. W. Drummond, 72 J. M. Izen, 72 X. C. Lou, 72 F. Bianchi afc , 73 D. Gamba ab , 73 M. Pelliccioni ab , 73 
M. Bomben afc , 74 L. Bosisio ab , 74 C. Cartaro ab , 74 G. Delia Ricca ab , 74 L. Lanceri ab , 74 L. Vitale ab , 74 V. Azzolini, 75 
N. Lopez-March, 75 F. Martinez- Vidal, 75 D. A. Milanes, 75 A. Oyanguren, 75 J. Albert, 76 Sw. Banerjee, 76 B. Bhuyan, 76 

H. H. F. Choi, 76 K. Hamano, 76 G. J. King, 76 R. Kowalewski, 76 M. J. Lewczuk, 76 I. M. Nugent, 76 J. M. Roney, 76 
R. J. Sobie, 76 T. J. Gershon, 77 P. F. Harrison, 77 J. Ilic, 77 T. E. Latham, 77 G. B. Mohanty, 77 E. M. T. Puccio, 77 

H. R. Band, 78 X. Chen, 78 S. Dasu, 78 K. T. Flood, 78 Y. Pan, 78 R. Prepost, 78 C. O. Vuosalo, 78 and S. L. Wu 78 

(The BABAR Collaboration) 

I Laboratoire d'Annecy-le-Vieux de Physique des Particules (LAPP), 
Universite de Savoie, CNRS/IN2P3, F-74941 Annecy-Le-Vieux, France 

2 Universitat de Barcelona, Facultat de Fisica, Departament ECM, E-08028 Barcelona, Spain 
:i INFN Sezione di Bari a ; Dipartimento di Fisica, Universitd di Bari b , 1-70126 Bari, Italy 
* University of Bergen, Institute of Physics, N-5007 Bergen, Norway 
5 Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 
6 University of Birmingham, Birmingham, B15 2TT, United Kingdom 
7 Ruhr Universitat Bochum, Institut fur Experimentalphysik 1, D-44780 Bochum, Germany 
8 University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1 
9 Brunei University, Uxbridge, Middlesex UB8 3PH, United Kingdom 
10 Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

II University of California at Irvine, Irvine, California 92697, USA 

12 University of California at Los Angeles, Los Angeles, California 90024, USA 
13 University of California at Riverside, Riverside, California 92521, USA 
lif University of California at San Diego, La Jolla, California 92093, USA 
15 University of California at Santa Barbara, Santa Barbara, California 93106, USA 
16 University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 
17 California Institute of Technology, Pasadena, California 91125, USA 

18 University of Cincinnati, Cincinnati, Ohio 45221, USA 

19 University of Colorado, Boulder, Colorado 80309, USA 

20 Colorado State University, Fort Collins, Colorado 80523, USA 
21 Technische Universitat Dortmund, Fakultdt Physik, D-44221 Dortmund, Germany 
22 Technische Universitat Dresden, Institut fur Kern- und Teilchenphysik, D-0 1062 Dresden, Germany 
23 Laboratoire Leprince-Ringuet, CNRS/IN2P3, Ecole Polytechnique, F-91128 Palaiseau, France 
■* University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 
25 INFN Sezione di Ferrara a ; Dipartimento di Fisica, Universitd di Ferrara b , 1-44100 Ferrara, Italy 

26 INFN Laboratori Nazionali di Frascati, 1-00044 Frascati, Italy 
27 INFN Sezione di Genova a ; Dipartimento di Fisica, Universitd di Genova b , 1-16146 Genova, Italy 

28 Harvard University, Cambridge, Massachusetts 02138, USA 
29 Universitat Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 
30 Humboldt- Universitat zu Berlin, Institut fur Physik, Newtonstr. 15, D- 12489 Berlin, Germany 
"imperial College London, London, SW7 2AZ, United Kingdom 
32 University of Iowa, Iowa City, Iowa 52242, USA 
33 Iowa State University, Ames, Iowa 50011-3160, USA 
Johns Hopkins University, Baltimore, Maryland 21218, USA 



3 



15 Laboratoire de lAccelerateur Lineaire, IN2P3/CNRS et Universite Paris-Sud 11, 
Centre Scientifique d'Orsay, B. P. 34, F-91898 Orsay Cedex, France 
36 Lawrence Livermore National Laboratory, Livermore, California 94550, USA 
37 University of Liverpool, Liverpool L69 7ZE, United Kingdom 
Queen Mary, University of London, London, El 4NS, United Kingdom 
39 University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 
40 University of Louisville, Louisville, Kentucky 40292, USA 
* Johannes Gutenberg- Universitat Mainz, Institut fur Kernphysik, D-55099 Mainz, Germany 
42 University of Manchester, Manchester M13 9PL, United Kingdom 

43 University of Maryland, College Park, Maryland 20742, USA 
44 University of Massachusetts, Amherst, Massachusetts 01003, USA 
45 Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, USA 
46 McGill University, Montreal, Quebec, Canada H3A 2T8 
47 INFN Sezione di Milano a ; Dipartimento di Fisica, Universita di Milano b , 1-20133 Milano, Italy 
48 University of Mississippi, University, Mississippi 38677, USA 
49 Universite de Montreal, Physique des Particules, Montreal, Quebec, Canada H3C 337 
50 Mount Holyoke College, South Hadley, Massachusetts 01075, USA 
51 INFN Sezione di Napoli a ; Dipartimento di Scienze Fisiche, 
Universita di Napoli Federico IF, 1-80126 Napoli, Italy 
52 NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands 

5 University of Notre Dame, Notre Dame, Indiana 46556, USA 

54 Ohio State University, Columbus, Ohio 43210, USA 

55 University of Oregon, Eugene, Oregon 97403, USA 

56 INFN Sezione di Padova a ; Dipartimento di Fisica, Universita di Padova b , 1-35131 Padova, Italy 
jl Laboratoire de Physique Nucleaire et de Hautes Energies, 
IN2P3/CNRS, Universite Pierre et Marie Curie-Paris6, 
Universite Denis Diderot- Paris7, F-75252 Paris, France 
58 University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 
59 INFN Sezione di Perugia a ; Dipartimento di Fisica, Universita di Perugia b , 1-06100 Perugia, Italy 
60 INFN Sezione di Pisa a ; Dipartimento di Fisica, 
Universita di Pisa b ; Scuola Normale Superiore di Pisa c , 1-56127 Pisa, Italy 
61 Princeton University, Princeton, New Jersey 08544, USA 
62 INFN Sezione di Roma a ; Dipartimento di Fisica, 
Universita di Roma La Sapienza b , 1-00185 Roma, Italy 
63 Universitat Rostock, D-18051 Rostock, Germany 
64 Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 
65 CEA, Irfu, SPP, Centre de Saclay, F-91191 Gif-sur-Yvette, France 
G6 SLAC National Accelerator Laboratory, Stanford, California 94309 USA 
67 University of South Carolina, Columbia, South Carolina 29208, USA 
68 Stanford University, Stanford, California 94305-4060, USA 
69 State University of New York, Albany, New York 12222, USA 

70 University of Tennessee, Knoxville, Tennessee 37996, USA 

71 University of Texas at Austin, Austin, Texas 78712, USA 
72 University of Texas at Dallas, Richardson, Texas 75083, USA 

7:1 INFN Sezione di Torino a ; Dipartimento di Fisica Sperimentale, Universita di Torino b , 1-10125 Torino, Italy 
74 INFN Sezione di Trieste 11 ; Dipartimento di Fisica, Universita di Trieste b , 1-34127 Trieste, Italy 
75 IFIC, Universitat de Valencia-CSIC, E-46071 Valencia, Spam 
76 University of Victoria, Victoria, British Columbia, Canada V8W 3P6 
77 Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 
78 University of Wisconsin, Madison, Wisconsin 53706, USA 
(Dated: June 2, 2009) 

We have performed a search for the rare leptonic decays B + — » £ + ue (I = e,fi), using data 
collected at the T(iS) resonance by the BABAR detector at the PEP-II storage ring. In a sample 
of 468 million BB pairs we find no evidence for a signal and set an upper limit on the branching 
fractions B(B + -» < 1-0 x iO -6 and B(B+ -> e+v e ) < 1.9 x 10" 6 at the 90% confidence level, 

using a Bayesian approach. 

PACS numbers: 13.20.-v, 13.25.Hw 

In the standard model (SM), the purely leptonic B SM branching fraction can be calculated as 0] 
meson decays B + — > l + vi [1[ proceed at lowest order 
through the annihilation diagram shown in Fig. [TJ The 
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B(B+ -» e+u t ) = 



GpmB'm'j 
8^ 



1 - -f /^iK 6 rr B , 

Win 



(1) 

where Gf is the Fermi coupling constant, mi and tub 
are respectively the lepton and B meson masses, and tb 
is the B + lifetime. The decay rate is sensitive to the 
CKM matrix element | V u b | [H and the B decay constant 
fs that describes the overlap of the quark wave functions 
within the meson. 




FIG. 1: Lowest order SM Feynman diagram for the purely 
leptonic decay B + — ■> l + vi. 

The SM estimate of the branching fraction for B + — > 
t+v t is (1.59 ±0.40) x 10~ 4 assuming t b = 1.638±0.011 
ps [3, V ub = (4.39±0.33)xl0~ 3 determined from in- 
clusive charmless semileptonic B decays [B| and Jb = 
216±22 MeV from lattice QCD calculation ||. To a very 
good approximation, helicity is conserved in B + — > fJ> Vn 
and B + — > e + v e decays, which are therefore suppressed 
by factors m 2 ^ e /m 2 with respect to B + — ► t + v t , lead- 
ing to expected branching fractions of B(B + — > /i + ^ A1 ) = 
(5.6±0.4)xlO- 7 and£(B+ -> e+z/ e ) = (1.3±0.4) x 10" 11 . 
However, reconstruction of B + — ► r + ^ r decays is ex- 
perimentally more challenging than B + —> y^v^ or 
B + — » e + i/ e due to the large missing momentum from 
multiple neutrinos in the final state. 

Purely leptonic B decays are sensitive to physics be- 
yond the SM, where additional heavy virtual parti- 
cles contribute to the annihilation processes. Charged 
Higgs boson effects may greatly enhance or suppress the 
branching fraction in some two-Higgs-doublet models Q • 
Similarly, there may be enhancements through mediation 
by leptoquarks in the Pati-Salam model of quark-lepton 
unification Direct tests of Yukawa interactions in 
and beyond the SM are possible in the study of these de- 
cays, as annihilation processes proceed through the lon- 
gitudinal component of the intermediate vector boson. 
In particular, in a SUSY scenario at large tan [3, non- 
standard effects in helicity-suppressed charged current 
interactions are potentially observable, being strongly 
tan /J-dependent and leading to Q : 



exp 



B{B- 



(1 - tan ; 



P M 2 H ] 



(2) 



Evidence for the first purely leptonic B decays has 
recently been presented by both the BABAR and Belle 
collaborations. The latest HFAG world average of the 
BABAR [| and Belle E3] results is B{B+ -> t+v t ) = 
(1.51 ±0.33) x 10~ 4 [11]. The current best published up- 
per limits on B + — > and B + — » e + v e are B(B + — » 
fi+v^) < 1.7 x 10~ 6 and B(B+ -» e + u e ) < 9.8 x 10~ 7 at 
90% confidence level from Belle using a data sample of 
235 fb" 1 [13]. 

The analysis described in herein is based on the en- 
tire dataset collected with the BABAR detector [l]| at 
the PEP-II storage ring at the T(AS) resonance ("on- 
resonance"), which consists of 468 million BB pairs, cor- 
responding to an integrated luminosity of 426 fb -1 . In 
order to study background from continuum events such 
as e + e~ — > qq (q = u,d,s,c) and e + e~ — > t + t~, an 
additional sample of about 41 fb^ 1 was collected at a 
center-of-mass (cm.) energy about 40 MeV below the 
T(4S) resonance ( "off- resonance" ) . 

In the BABAR detector, charged particle trajectories 
are measured with a 5-layer double-sided silicon vertex 
tracker and a 40-layer drift chamber, which are con- 
tained in the 1.5 T magnetic field of a superconducting 
solenoid. A detector of internally reflected Cherenkov 
radiation provides identification of charged kaons and pi- 
ons. The energies and trajectories of neutral particles are 
measured by an electromagnetic calorimeter consisting of 
6580 CsI(Tl) crystals. The flux return of the solenoid 
is instrumented with resistive plate chambers and, more 
recently, limited streamer tubes in order to provide 
muon identification. A GEANT4-based [HI Monte Carlo 
(MC) simulation of generic BB, qq, q — u,d, s, c, and 
t + t~ events as well as B + — > /i+i^ and B + — > e + v e 
signal events is used to model the detector response and 
test the analysis technique. 

The B + — > l + V£ decay produces a mono-energetic 
charged lepton in the B rest frame with a momentum 
p* w ms/2. The B mesons produced in T(45) decays 
have a cm. momentum of about 320 MeV/c, so we ini- 
tially select lepton candidates with cm. momentum 2.4 
< Pc.m. < 3.2 GeV/c, to take into account the smearing 
due to the motion of the B. A tight particle identification 
requirement is applied to the candidate lepton in order 
to discard fake muons or electrons. 

Since the neutrino produced in the signal decay is not 
detected, all charged tracks besides the signal lepton and 
all neutral energy deposits in the calorimeter are com- 
bined to reconstruct the companion (tag) B. We in- 
clude all neutral calorimeter clusters with cluster energy 
greater than 30 MeV. Particle identification is applied 
to the charged tracks to identify electrons, muons, pions, 
kaons and protons in order to assign the most likely mass 
hypothesis to each B tag daughter and thus improve the 
reconstruction of the B tag . Events which have additional 
lepton candidates are discarded. These typically arise 
from semileptonic B tag or charm decays and indicate the 
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presence of additional neutrinos, for which the inclusive 
Stag reconstruction is not expected to work well. 

The signal lepton's momentum in the signal B rest 
frame p* is refined using the -B tag momentum direction. 
We assume that the signal B has a cm. momentum of 
320 MeV/c and choose its direction as opposite that of 
the reconstructed B tag to boost the lepton candidate into 
the signal B rest frame. 

Signal events are selected using the kinematic variables 
AE = Eb — -Ebcam i where Eb is the energy of the -B tag 
and -©beam is the beam energy, all in the cm. frame. 
For signal events in which all decay products of the Stag 
are reconstructed, we expect the AE distribution to peak 
near zero. However, we are often unable to reconstruct all 
Stag decay products, which biases the AE distribution 
toward negative values. For continuum backgrounds, AE 
is shifted toward relatively large positive values since too 
much energy is attributed to the nominal Stag decay, 
while there is a negative bias in t + t~ events due to the 
unreconstructed neutrinos. 

We require the tag B to satisfy -2.25 < AE < GeV 
for B + — » [i^v^ decays. For B + — > e + v e decays, we 
require a linear combination of AE and the tag B trans- 
verse momentum pr to satisfy (pt + 0.529 • AE) <0.2 
and (pt — 0.529 • AE) <1.5. This selection rejects back- 
ground events arising from two-photon process e + e~ — > 
e + e~7*7*, 7*7* — > hadrons with one of the final elec- 
trons scattered at a large angle and detected. The coef- 
ficient of the AE term is extracted from the data. 

Backgrounds may arise from any process producing 
charged tracks in the momentum range of the signal, par- 
ticularly if the charged tracks are leptons. The two most 
significant backgrounds are B semileptonic decays involv- 
ing b — > ulvi transitions in which the momentum of the 
leptons at the endpoint of the spectrum approaches that 
of the signal, and from continuum and t + t~ events in 
which a charged pion is mistakenly identified as a muon 
or an electron. 

Continuum events tend to produce a jet-like event 
topology, while BB events tend to be more isotropically 
distributed in the cm. frame, and are suppressed us- 
ing event shape parameters. Five different spatial and 
kinematical variables, considered separately for B + — » 
p vv^ and B + — > e + i/ ei are combined in Fisher discrimi- 
nants [la ] . The most effective discriminating parameters 
are the ratio of the second L2 and the zeroth Lq mono- 
mial L n = Si\pi\ cos(a) ra , where the sum runs over all 
Stag daughters having momenta pt and a is the angle 
with respect to the lepton candidate momentum, both in 

the cm. frame, and the sphericity S = |min ^77^35 , 
where the T subscript denotes the momentum compo- 
nent transverse to the sphericity axis, which is the axis 
that minimizes S. S, in fact, tends to be closer to 1 
for spherical events and for jet-like events. In order to 
take into account the changes in detector performance 




FIG. 2: Distributions of signal (a,b) and background (c,d) 
rriES (left) and pfit (right) for B + — > fi + Vfi from MC simu- 
lation (a,b and c) and from tues sideband 5.17 < rriES < 5.2 
GeV/c 2 (d). 




I, (GeV/c 2 ) P FIT (GeV/c) 

FIG. 3: Distributions of signal (a,b) and background (c,d) 
rriES (left) and pfit (right) for B + — ■> e + u e from MC simula- 
tion. 



throughout the years, in particular in muon identifica- 
tion, the data sample is divided into six different data 
taking periods, and the Fisher discriminants and selec- 
tion criteria are optimized separately with the algorithm 
described in [ItJ for each period. 

The two-body kinematics of the signal decay is ex- 
ploited by combining the signal lepton momentum in the 
B rest frame p* and p c .m. in a second Fisher discrim- 
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inant (pfit) which discriminates against the remaining 
semileptonic bb and continuum background events which 
populate the end of the lepton spectrum in both frames. 
The p* and p c . m , coefficients in the linear combination are 



determined separately for B + 
with 



13- 



and B 



We employ an extended maximum likelihood (ML) 
fit to extract signal and background yields using si- 
multaneously the distributions of the Fisher output 
Pfit and the energy-substituted mass mESj defined as 



beam 



\ps I 2 , where %>b is the momentum of the re- 



constructed Stag candidate in the cm. frame. 

Signal hues and Pfyt probability density functions 
(PDFs) are fixed in the final fit and are parameterized 
from simulated events, respectively, with a Crystal Ball 
function [l8j and the sum of two Gaussians (double Gaus- 
sian) for both B + — > and B + — > e + v e . 

The background tties distribution is described by an 
ARGUS function whose slope is determined in the fit 
to the yields [19]. To parameterize the background 
Pfit distributions, we studied the possibility of using 
the niES sideband of on-resonance data. We found the 
B + — > /i + z^ sideband suited for this purpose, while the 
B + — > e + v e sideband is not sufficiently populated. We 
use the region 5.17 < tues < 5.2 GeV/c 2 to parameter- 



ize the B + 



background pfit distribution, and 



simulated events for the background B + — > e + v e pfit 
distribution. 

Separately for B + — > fJ^v^ and B + — > e + v e , the sum 
of two Gaussians with different sigmas on the right and 
the left of the mean (bifurcated Gaussians) is used to 
parameterize the background pyyt distribution and the 
relative fraction of the two bifurcated Gaussians is de- 
termined from the fit to the data. Figures [2] and [3] show 
background and signal tues and pfit distributions for 
B + — > [i^v^ and B + — > e + z/ e , respectively, with the 
PDFs described above superimposed. 

In the on-resonance data the ML fit returns 1 ± 15 
signal B + — > fi^v^ candidate events and 18 ± 14 signal 
B + — > e + v e candidate events. Distributions of the fit 
data events with the final fit superimposed, as well as 
the signal and background PDFs, are shown in Figured] 
for B + — > ix + v^ and B + — > e + v e , respectively, projected 
on ities and pfit- 

We next evaluate systematic uncertainties on the num- 
ber of B^ in the sample, the signal efficiency and the sig- 
nal yield. The number of B^ mesons in the on-resonance 
data sample is estimated to be 468 x 10 6 with an uncer- 
tainty of 1.1% [20| , as suming equal B + and B° produc- 
tion at the T(45) (2l|. 

The uncertainty in the signal efficiency includes the 
lepton candidate selection (particle identification, track- 
ing efficiency and event selection Fisher requirement) as 
well as the reconstruction efficiency of the tag B. The 
systematic uncertainty on the particle identification ef- 




(b) 



(.2 5.22 5.24 5.26 5.28 5.3 
nijjs [GeV/c 2 ] 



10 15 20 

p (GeV/c) 




ficiency is evaluated using e + e 



FIG. 4: Final fit to the data projected on ities (left) and 
Pfit (right) distributions for B + — > /-L + u^ events (a,b) and 
B + — > e + u e events (c,d) : the solid blue line is the total PDF, 
the dashed red line is the background PDF and the dashed- 
dotted black line is the signal PDF 



e + e~[i + /i~ and Bhabha event control samples derived 
from the data, which are weighted to reproduce the kine- 
matic distribution of the lepton signal candidate. Com- 
paring the cumulative signal efficiency obtained with and 
without these weights, a total discrepancy of 1.9% for 
B + — » and 2.3% for B + — > e + u e is found and this 

value is taken as the particle identification systematic 
uncertainty. Tracking efficiency is studied employing t 
decays, which must produce an odd number of final state 
charged tracks because of charge conservation. Thus, one 
can determine an absolute efficiency because the number 
of events with a missing track can be measured. The 
uncertainty associated with the tracking efficiency and 
the data/MC discrepancy evaluated with this method are 
taken in quadrature for a total tracking efficiency uncer- 
tainty of 0.4% per track. 

In order to evaluate the systematic uncertainty as- 
sociated with the requirements on the Fisher discrimi- 
nants, we compare data and MC Fisher distributions in 
the sidebands AElO for the B + — > n + v^ sample and 
(p T + 0.529 • AE)i 0.2 for the B+ -> e+v e sample. We 
fit the data/MC ratio with a linear function, with re- 
sults consistent with a unitary ratio in the whole Fishers 
range. We take the error on the intercept as the system- 
atic uncertainty on the Fisher discriminants, that is 1.4 
% for B+ -> n+Vp and 5.3% for B+ -> e+v e . 

The tag B reconstruction has been studied with a con- 
trol sample of B + — > Z?(*)°7r + events, where the D is 
reconstructed into D° —> K + ix~ and D° — > K~it + , and 
the D* into D*° -> D°j or D*° -> L>°7r°. These two- 



7 



TABLE I: Contributions to the systematic uncertainty on 
the signal efficiency. Total systematic represent the sum in 
quadrature of the table entries. 



Source 




B+ - 


e + v e 


Particle identification 


1.9% 


2.3 


% 


Tracking efficiency 


0.4% 


0.4 


% 


Tag B reconstruction 


3.0% 


0.4 


% 


Fisher selection 


1.4% 


5.3 


% 


Total 


3.8% 


5.8 


% 



body decays are topologically very similar to our signal, 
as the charged pion can be treated as the signal lepton 
and the £)(*)° decays products ignored to simulate the 
missing neutrino. The tag B reconstructed in the con- 
trol sample thus simulates the tag B reconstruction in the 
nominal data sample. We compare the efficiencies for our 
tag B selection cuts in the B+ -> £>M°7r + data and MC 
to quantify any data/MC disagreements that may affect 
the signal efficiency. We find a data/MC discrepancy on 
B+ -> D(*)°7T+ control sample of 3.0% for B+ -> /i+z^ 
decays and 0.4% for B + — > e + v e decays, and assign these 
as the signal efficiency uncertainty arising from the tag 
B selection. 

A summary of the systematic uncertainties in the sig- 
nal efficiency is given in Table Q] The final B + — > fi + v^ 
signal efficiency is (6.1 ± 0.2)% and the B + — > e + v e sig- 
nal efficiency is (4.7 ± 0.3)%, where the errors are the 
sum in quadrature of statistical and systematic uncer- 
tainties. 

The systematic uncertainty in the yields comes from 
the pfit and rriES PDF parameters, which are kept fixed 
in the final fit and, in the B + — > e + v e case, from the 
use of MC simulation to extract the PDF shapes. The 
fit parameters extracted from MC are affected by an 
uncertainty due to MC statistics. In order to evaluate 
the systematic uncertainty associated with the parame- 
terization, the final fit has been repeated 500 times for 
each background and signal PDF parameter which is kept 
fixed in the final fit. We randomly generate the PDF pa- 
rameters assuming Gaussian errors and taking into ac- 
count all the correlations between them. We perform a 
Gaussian fit to the distribution of the number of signal 
events for each parameter, take the fitted sigma as the 
systematic uncertainty and sum in quadrature. The to- 
tal systematic uncertainty on the signal yield from all 
signal and background PDF parameters is 8 events for 
B + — » and 10 events for B + — > e + v e . 

For the B + —> e + v e sample, an additional systematic 
uncertainty coming from possible discrepancies in the 
shape of the pfit background distribution in data and 
simulated events must be accounted for. The data/MC 
ratio of the pfit distribution in the rriES sideband 5.16 
< rriES < 5.22 GeV/c 2 is fit with a linear function. The 



background pfit distribution shape is varied according 
to the fitted linear function and its associated statistical 
uncertainties; the total systematic contribution from this 
procedure is 4 events. 

To evaluate the branching fraction we use the following 
expression: 

B{B^l+v) UL = -^L_ (3) 
N B ± ■ e 

where N sig represents the observed signal yield, N B ± the 
number of B + B~ in the sample (where equal production 
of B + B~ and B°B a is assumed) and e is the signal effi- 
ciency. 

As we did not find evidence for signal events, we em- 
ploy a Bayesian approach to set upper limits on the 
branching fractions. Flat prior in the branching frac- 
tions are assumed for positive values of the branching 
fractions and Gaussian likelihoods are adopted for the ob- 
served signal yield, related to B by Eq.([3]). The Gaussian 
widths are fixed to the sum in quadrature of the statisti- 
cal and systematic yield errors. The effect of systematic 
uncertainties associated with the efficiencies, modeled by 
Gaussian PDFs, is taken into account as well. We extract 
the following 90 % confidence level upper limits on the 
branching fractions: 

B(B+ -> (j+v,,.) < 1.0 x 10~ 6 (4) 
B(B + -> e + v e ) < 1.9 x 10~ 6 . (5) 

The 95% upper limits are B{B+ -> /i+i^) < 1.3 x 10~ 6 
and B(B + — > e + v e ) < 2.2 x 10 -6 . This result improves 
the previous best published limit for B + — > branch- 
ing fraction by nearly a factor two, to a value twice the 
SM prediction. The B + — > e + v e result is consistent with 
previous measurements. It should be noted that the re- 
sults in are obtained using a different statistical ap- 
proach to interpret the observed number of signal events. 
The results show no deviation from the SM expectations. 
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